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Abstract 

We study the tunneling magneto-transport properties of the Ferromagnetic Insulator-Normal Insulator-Ferromagnetic 
Insulator(F|N|F) and Ferromagnetic Insulator-Barrier Insulator-Ferromagnetic Insulator (F|B|F) junctions on the surface of topo- 
logical insulator in which in-plane magnetization directions of both ferromagnetic sides can be parallel and antiparallel. We derive 
analytical expressions for electronic conductances of the two mentioned junctions with both parallel and antiparallel directions of 
magnetization and using them calculate magnetoresistance of the two junctions. We use thin barrier approximation for investigating 
the F|B|F junction. We find that although magnetoresistance of the F|N|F and F|B|F junctions are tunable by changing the strength 
'"^ of magnetization texture, they show different behaviors with variation of magnetization. In contrast to the magnetoresistance of 
F|N|F, magnetoresistance of F|B|F junctions shows very smooth enhance by increasing the strength of magnetization. We suggest 
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an experimental set up to detect our predicted effects. 
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Q_il. Introduction 

^ Recently, some topological insulators (TI) are observed ex- 
"£3 perimentally and investigated theoretically fl] [2] [3] @] |7] 
g El |3 El HI [TTJ [U. A TI is an insulator in bulk but a ro- 
i bust conducting in edge even in presence of deformation or 
*^ strong disorder i.e. it is gapless only in surface. So it can 
^ carry charge current on its surface even in the absence of mag- 
^ netic field. This state of condensed matter is due to spin-orbit 
i— —i coupling and protected by origin of time-reversal symmetry 
019] 0. For two-dimensional TI we have helical modes (the 
►^J spin orientation of surface electrons are perpendicular to mo- 
mentum direction ) at the edges [2|, similar to quantum Hall 
00 systems. This property is reason for some interesting phys- 
00 ical properties of the system|10|. For three-dimensional TI, 
N J # strong (weak) TI is defined for odd (even) number of chiral 
Dirac fermions on surface 0. In graphene, There are four 
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Dirac cones in momentum space [ 15], but the number of Dirac 
cones in TI's edge are odd. So there are good unconventional 
properties such as fractional quantum Hall effect, n Berry's 
phases and immunity to Anderson localization due to odd spin- 
Dirac on the surface ifTl l6l [161 flTl . On the other hand, study- 
ing of transport properties is useful to understand the macro- 
scopic and applicatory properties. Some attractive physical 
properties are seen when TI is on proximity of ferromagnet or 
superconductor lfTSl [19l l20l [131 [T4l . Also the tunneling magne- 
toresistance (TMR) effect has attracted much interest because 
of crucial applications in spintronics ETl . In this paper we in- 
vestigate TMR and transport properties of the Dirac fermions 
on the surface of topological insulator for F|N|F and F|B|F junc- 
tions within the clean limit and low temperature regime. We 
assume that the used topological insulator is made of HgTe or 




^ 1 Topological Insulator 



Figure 1: (Color online) Schematic model of the suggested set up for a) F|N|F 
and h) F|B|F junctions on top of the topological insulator. We assume that mid- 
dle of the junctions is located at x = d/2 and directions of the magnetizations 
in both sides can be parallel and antiparallel. For F|B|F case we employe thin 
barrier approximation namely, when transverse voltage Vo — > oo and d — > oo. 
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Figure 2: (Color online) The electronic conductances for parallel (a) and an- 
tiparallel (b) directions of magnetization in F|N|F junctions. Here we have fixed 
the M/Mc = 0.9 and plotted electronic conductances vs. width of the normal 
region (KpL = ^-L) for three values of eV gate voltage. 



Bi2Sej. The advantage of the mentioned alloys is their single 
Dirac cone in Berillouin zone in which two dimensional Dirac 
equation govern the Dirac fermions on surface of the topologi- 
cal insulator|[T]|2]|6]|- Due to the exchange coupling, the ferro- 
magnetism can be induced on the surface of topological insula- 
tor when a ferromagnetic layer is deposited on top of it[24|. We 
use two dimensional Dirac Hamiltonian and derive wave func- 
tions in all regions. Using boundary conditions we obtain ana- 
lytical expressions for reflection and transmission probabilities. 
Using the obtained reflection and transmission probabilities, we 
calculate electronic conductance of the two mentioned systems 
when the magnetization of ferromagnetic layers are parallel and 
antiparallel. To investigate properties of F|B|F junctions we use 
the thin barrier approximation. The used barrier can be made by 
applying a transverse voltage or doping 1 22 23 ]. The main con- 
cern is to investigate the magnetoresistance of the mentioned 
systems and their differences. 
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Figure 3: (Color online) The magnetoresistance of the F|N|F junctions vs. width 
of the normal region (KpL = -^rpL) for three values of M/Mc. eV the gate 
voltage of junction is fixed at 0.1/j. 



Figure 4: (Color online) The magnetoresistance of the F|N|F junctions vs. 
M/M c for three values of gate voltage eV and the width of the junction is fixed 
at KpL=10, 



2. Tunneling Magnetoresistance in F|N|F Junctions 

On the surface of topological insulator Dirac fermions are 
governed by two dimensional Dirac Hamiltonian that is given 
by: 



H = hvp&x — fil 



(1) 



where & is vector of Pauli matrices in spin space and I is iden- 
tity matrix. Also, vf is Fermi velocity and /i is chemical po- 
tential. When a ferromagnetic electrode with magnetization m 
is deposited on surface of the TT, due to the proximity of TI 
and ferromagnetism, an additional term appear into Eq.([T]i as 
follow! 24]: 



Hq = hvpCT.m. 



(2) 



We consider an F|N|F junction in which thickness of normal re- 
gion is d. We assume the TI is located in xy-plane and two 
ferromagnetic electrodes is deposited on TI such that m = 
ymo{9(-x) + 9(x - d)) in which +(-) refers to parallel (antiparal- 
lel) in magnetization alignment of electrodes (See Fig. [TJ. Here 
0(x) is the well known step function. By solving the Eq. ([T| for 
x < region, Dirac spinors read as: 
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(±ik x x+ik v y) 
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±e ±ip 



(3) 



where k x and k y are components of wave vector in x and y di- 
rections, respectively. Also the + sings refer to the directions of 
the x-axis and /3 is incidental angle of the Dirac fermions which 
j6 and k x are described as: 

%v F {k y + M) 



arcsin 



\e + fi\ 



^ 2 
nvp 



(k y + M) 2 , 



(4) 



where M — mo/fiVf depends on the exchange coupling and one 
can tune the parameter via an external field[24]. For normal 
region between < x < d, Dirac spinors are given by: 



^n(x,y) = —e' 



(±ik" x+ik v y) 



1 

±e ±: 



(5) 



in which k" is x-component of the wave vector and a is angle 
of incidental Dirac fermions and are defined as: 
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In region x > d, there are two choices, parallel or antiparallel 
directions for magnetization of the electrodes. When we make 
our system with parallel magnetization, all spinors are similar 
to Eq. <[3j but for antiparallel case the spinors read: 



o < 



(b) 
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l 



(±ik[x+ik y y) 
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in which B' and k' x are defined as: 

TM^y - M) 



Figure 5: (Color online) The electronic conductances of F|B|F junctions within 
the thin barrier approximation vs. strength of barrier^ = Vod/TiVf for three 
values of gate voltage eV. Here M/Mc is fixed at 0.9. 
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As seen in Eqs.([4} [8j, for antiparallel magnetization case when 
M is larger than its critical value namely Mq = \e + /j]/Hvf, the 
wave function in Eq. (|7]i decays. For an incidental particle from 
x < into the F|N interface at x = 0, the transverse wave vector 
k y and excitation energy e are conserved in scattering process. 
The Dirac spinors within each region must satisfy the following 
boundary conditions at x — 0: 



xf/ 1 + nffj = a\\ft n + a 2 ifr u . 
and at x = d : 

t\\i\ = anffu + a 2 tfr n . 



(9) 



(10) 
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By solving Eqsj9]and 10 transmission probability T = \t\ 2 for 
parallel and antiparallel magnetization directions are derived as: 



Figure 6: (Color online) The magnetoresistance of the F|B|F junctions within 
the thin approximation vs. strength of barrier^' = Vgd/tiVf for three values of 
M/Mc and fixed gate voltage eV = 0.1/j. 
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cos fi cos a 



A cosH^f) cos 2 a + B [cos(^) - sin a sin(^)] 2 



(11) 



cos B cos a 



A cos 2 (B) cos 2 a + B [1 - sinasin(jS)] 2 



(12) 



where A = cos 2 (-^M) and B = sin z (^M)- The P and AP 
indexes stand for parallel and antiparallel magnetizations. The 
EqsfTT|and 12 are the main results of this section. Using the 
transmission probabilities and conductance formula[23], one 
can calculate the conductances of the mentioned structures con- 
taining parallel and antiparallel magnetizations, 



G(e) = ^ f ' TcosBdB 

1 J-n/2 



(13) 



here Go = N(e)We 2 /nft 2 Vf is normal conductance and A^(e) = 
|/i + eV\/2n(hvf ) 2 is density of state of Dirac fermions and W is 
width of the junction. 

The conductances of the parallel and antiparallel magnetization 
are plotted versus kfL in Fig|2j here kf = fi/hvf is Fermi vector 



of the system. The oscillatory behavior of both conductances 
are results of coherent interference of the Dirac fermions at two 
interfaces. By using the conductances of parallel and antipar- 
allel magnetization, one can calculate the magnetoresistance of 
the system which is defined as MR = 100% x (Gp - Gap)IGp. 
This quantity may have many interesting and important appli- 
cations in spintronic industry, in similarity with the Ref. |21 1. 
The plotted magnetoresistance as a function of kfL is shown in 
Fig(3] which shows oscillatory behavior for all of the gate volt- 
ages as we expected. This oscillatory behavior is a result of 
oscillatory behavior in conductances for parallel and antiparal- 
lel magnetizations. The plotted magnetoresistance as a function 
of M can be seen in Fig|4]that is also one of the main results of 
this section. By increasing M, the magnetoresistance enhance 
to reach the perfect value in Mq- Although for small values of 
M the magnetoresistance oscillates, the quantity for large val- 
ues of M, shows a smooth behavior. Now we proceed to present 
an investigation of the FI/BI/FI junction. 



paper. As seen in Eq. ( 18 1, for an incidental particle with nor- 
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Figure 7: (Color online) The magnetoresistance of the F|B|F junctions within 
the thin approximation vs. M/Mc for three values of strength of barrier^ = 
Vod/Hvf and fixed gate voltage eV = 0.1//. 



3. Tunneling Magnetoresistance in F|B|F Junctions 

Now we focus our attention to F|B|F junctions. For pro- 
ducing such structures one can apply a transverse gate volt- 
age into the normal region and make it as a barrier for moving 
Dirac fermions from left ferromagnet electrode to interface (See 
Fig[T[ part (b)). By applying the barrier voltage, another addi- 
tional term added to the Dirac Hamiltonian Eq. ([T} for normal 
region and change it as: 



H = hv F &X + Vol - ill. 
Dirac fermion states within the barrier region are: 
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(14) 
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where 6 and k" are defined as: 
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We are interested in the influence of thin barrier regime on the 
magnetoresistance. In this limit, width of the junction d ap- 
proaches to zero and barrier voltage Vq move to large values 
simultaneously so that the barrier strength x — eVod/hvf re- 
mains constant. In this limit 8^0 and k'^ d ~ X- Substituting 



the present Dirac spinors Eq. (15 1 into the boundary condi 



tions Eqs|9]and 10 , we have derived the transmission probabil- 
ities for parallel (P label) and antiparallel (AP label) structures 
within the thin barrier regime as below: 



' AP 



COS B 



A' cos 2 (^f) + B' [cos(£/)-sin(^)] 2 
cos 2 /? 

A' cos 2 (J3) + B' [1 - sin(/3)] 2 ' 



(17) 



(18) 



where A' 



cos 2 (-^) and B' 



sin (^). The above obtained 



analytical expressions for the transition probabilities of paral- 
lel and antiparallel structures are other results of the present 



mal direction to the interface, the transmission probability is 
equal to unity. This situation is similar to the Klein tunneling in 
quantum relativistic theory ||25l . To calculate the conductance 
of the mentioned parallel and antiparallel structure within the 
thin barrier regime, one should use the obtained transmission 



probabilities (Eqs 17|18 i in Eq. ( 13 1. The plotted conductances 
for both parallel and antiparallel magnetization are shown in 
Fig|5]as a function of the barrier strength^. Changing the gate 
voltage doesn't make any changes in the positions of the ex- 
ternum in conductances of parallel and antiparallel magneti- 
zation. The maximum of parallel conductance has n/2 phase 
difference with respect to antiparallel case. Unlike damping 
oscillatory behavior of MR versus length of the normal region 
KfL in F|N|F case (See Fig(3]), here the oscillatory behavior of 
the conductance leads to harmonic oscillations in MR, as shown 
in Figj6] We plotted the magnetoresistance for this set up as a 
function of M, as shown in FigJ7| which approximately shows 
linear behavior. This linear manner in magnetoresistance for 
F|B|F junction seems to be more applicatory than the discussed 
behavior of magnetoresistance in F|N|F junctions. By tuning 
the gate voltage and exchange coupling due to the ferromag- 
netic electrodes, we can access to an arbitrary value in magne- 
toresistance which is useful in spintronics. 

4. Conclusions 

In summary, we have investigated the electronic trans- 
port properties of topological insulator-based ferromag- 
netic/normal/ferromagnetic (F|N|F) and topological insulator- 
based ferromagnetic/barrier/ferromagnetic F|B|F junctions on 
the surface of topological insulator in low temperature and 
clean regimes. We assumed the two magnetizations are in-plane 
and deposited on the surface of topological insulator. We have 
derived analytical expressions for transmission probabilities of 
two different magnetization directions, parallel and antiparallel 
for F|N|F and F|B|F junctions. To investigate electronic proper- 
ties of the F|B|F junction, in our calculations we assumed that 
d = and Vq — > co simultaneously. This approximation is well 
known as thin barrier approximation. Using the derived transi- 
tion probabilities, conductances and magnetoresistances of the 
systems were calculated. We have found that the very important 
quantity of such junctions (magnetoresistance) shows different 
behavior by increasing the magnetization strength of ferromag- 
netic region for F|N|F and F|B|F junctions. Although magne- 
toresistance of the two junctions can be tuned from values near 
up to 100%, magnetoresistance for F|N|F junctions shows os- 
cillatory behavior vs. M/Mc the magnetization strength, varia- 
tion of the quantity vs. M/M c for F|B|F junctions behave like a 
very smooth enhancement. 
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